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ABSTRACT 

We present the internal kinematics of UCD3, the brightest known ultra-compact dwarf galaxy 
(UCD) in the Fornax cluster, making this the first UCD with spatially resolved spectroscopy. 
Our study is based on seeing-limited observations obtained with the ARGUS Integral Field 
Unit of the VLT/FLAMES spectrograph under excellent seeing conditions (Of' 5 - Of' 67 
FWHM). 

The velocity field of UCD3 shows the signature of weak rotation, comparable to that 
found in massive globular clusters. Its velocity dispersion profile is fully consistent with an 
isotropic velocity distribution and the assumption that mass follows the light distribution 
obtained from Hubble Space Telescope imaging. In particular, there is no evidence for the 
presence of an extended dark matter halo contributing a significant (^33 per cent within 
R < 200 pc) mass fraction, nor for a central black hole more massive than ~5 per cent of the 
UCD's mass. While this result does not exclude a galaxian origin for UCD3, we conclude that 
its internal kinematics are fully consistent with it being a massive star cluster. 

Key words: galaxies: dwarf - galaxies: kinematics and dynamics - galaxies: evolution - 
galaxies: star clusters 



1 INTRODUCTION 

During the past decade, the existence of a previously unknown 
morphological class of compact stellar systems, called ultra- 
compact dwarf galaxies (UCDs), has been established. Objects of 
this type were first discovered in the dense environments of the 
nearby galaxy clusters like Fornax (Hilker et al. 1999, D rinkwa-"] 
|ter et al.||2000| ) and Virgo ( |Ha§egan et al.||2 005). They also were 
found around the dominant giant galaxies in nearby groups, e.g. 
NGC 5128 ( [Rejkuba et al.|2007l|Taylor et al.|2010|) and M 104 ([Haul 
|et al.|200 9^, as well as in compact groups (Da R ocha et al.|201l| K 

UCDs are characterised by predominantly old stellar popula- 
tions (> 8 Gyr, e.g. [Chilingarian et al.|2011| ), typical luminosities 
of - 14.0<My<- 10.5 mag, half-light radii of 10 < r/, < 100 pc 
(e.g. [Evstigneeva et al.||2007| ) and central velocity dispersions of 
20 < Oq < 50 km s~ ^ , yielding dynamical masses of 2 x 10^ <m< 
10^ M0 (e.g. Hilker et al. 2007 , Mieske et al. 2008 ) . Hence, they 
are larger, brighter and more massive than the most massive Milky 
Way globular clusters (GCs). At the same time they are signifi- 
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cantly more compact than dwarf galaxies of comparable luminosity 
( [Dabringhausen et al.|2008| ). 

One of the most striking ensemble properties of UCDs is that 
on average, their dynamical mass to light ratios (M/L) are ^30- 
50 per cent higher than the predictions of canonical stellar pop- 
ulation models (e.g. [Kruijssen & Mieske||2009| ). Moreover, their 
V-band M/Ly are on average twice as large as those of Galactic 
GCs of comparable metallicity (e.g. Mieske et al. 2008 , jDabring-| 
[hausen et"a l. 2008, Taylor et al. 2010). Depending on their origin, 
this may indicate that UCDs mark the onset of dark matter domi- 
nation in small stellar systems ( [Baumgardt & Mieske|2008 JGoerdt 



et al. 



2009 



2008), or harbour massive central black holes fMerritt et al. 
or exhibit an unusual initial stellar mass function (Mieske & 



Kroupa|2008l|Dabringhausen et al.|2009|[20T0l ) 



The nature of UCDs is still uncertain. Do they represent the 
most massive star clusters (Fellhauer & Kroupa 2 002 , Mieske et al.| 
2004 )? Or are they the result of the environmental transformation 
of galaxies ( |Bekki et al.|2003] )? 

Spatially resolved kinematics of UCDs open up a way to new 
insights on this matter, complementary to those gained from the 
study of the spatial and velocity distributions of UCD populations 
as a whole (e.g. Mieske et al.||2004|), o f structural properties ( |De| 
[Propris et al.||2005| [Evstigneeva et al. "2008), of integrated inter- 
nal velocity dispersions (e.g. Hilker et al._2007j , or of their stellar 
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Figure 1. Left: The upper row shows the ACS/HRC F606W image, the image reconstructed from the merged data cube with an overlaid polar grid used for the 
spatial binning and the weight-map, which illustrates the location of bad spaxels in the input data cubes. The bottom row shows the mean signal to noise ratio 
in each spatial bin per spectral bin (of 1 1 km s~^ or ~ 0.2 A), the recovered velocity field with a line indicating the axis of rotation and the velocity dispersion 
map. The scale of the images is indicated by the 2" scale bar at the lower left, \" corresponds to ~ 92 pc at the distance of the Fornax cluster. Right: Two 
example spectra (black curves), the best-fitting spectral models (blue), residuals (grey) and noise models (green lines). The red line represents the third order 
multiplicative polynomial used in the fitting to absorb the continuum normalisation. The two spectra are representative for the range of quality in our data, 
corresponding to a central spatial bin (top) and an outer spatial bin, which is affected by the calibration lamp stray light (bottom). 



population content (e.g. l Taylor et al.|2010[[Paudel et al.|20T0l ). Re- 
solving the kinematics of objects with angular half-light radii on 
the order of I" , however, takes optical integral field spectroscopy 
to its limits. On the other hand, present day adaptive-optics aided, 
near infra-red integral field units (IFUs; e.g. Gemini/NIFS, |Mc-| 
IGregor et al.|[2003l VLT/SINFONI, |Eisenhauer "erari |2003] ), pro- 
vide only moderate spectral resolution (R = AA/A < 4000 or 
<^ins ^ 32 km s~^), limiting the ability to accurately sample the 
velocity distribution of typical UCDs with velocity dispersions of 
~ 25 km s~^. 

In this Letter we present, for the first time, the spatially re- 
solved kinematics of an UCD. A detailed description of the data 
reduction, of the technical aspects involved in the analysis and of 
the dynamical modelling will be presented in a forthcoming paper. 
The target of our study, UCDS, is one of the originally discovered 
objects (Hilk er et al.|19 99 ), which subsequently led to the defini- 
tion of the class of UCDs. Nevertheless, UCDS is not an average 
UCD in several ways: With My = —IS. 55 mag, it is the brightest 
known UCD in the Fornax cluster, separated by ^ 0.8 mag from the 
rest of the cluster's UCD population. In projection, it is relatively 
close to the central cluster galaxies NGC 1S99 and NGC 1404 
(projected distances of 46kpc and ISkpc, respectively). It is un- 
usually red, with a spectroscopic metallicity of [Fe/H]~ — 0.2dex 
( [Chilingarian et al. ] |2008| ). Most importantly, it has a large half- 
light radius of = 87 pc (or ^ 0'/95) and is one of the few UCDs 
where an extended, low surface-brightness envelope has been de- 
tected ( [Evstigneeva et al.||2008| ), making it possible to resolve it 
from the ground with seeing-limited instruments. 



2 OBSERVATIONS AND DATA REDUCTION 

UCDS was observed under excellent seeing with the ARGUS IFU 
( |Kaufer et al.[200Sj of the Flames/GIRAFFE ( |Pasquini et al.|2002| ) 



spectrograph mounted on UT2 at the VLT. The IFU was used in 
the 1:1 scale, yielding a field of view of llf'5 x 7f'S at a spatial 
sampling of 0^/52 x 0^/52 per lenslet (or spaxel), with the LR04 
grism, providing a spectral resolution of R~9600 (or ^ IS km s~^ 
in terms of Gaussian a) over the range 5015 — 58S1 A. The obser- 
vations were executed in service mode over the period of almost 
one year and consisted of eight exposures of 2775 seconds each. 

The basic data reduction and the extraction of wavelength cal- 
ibrated spectra was carried out using the Giraffe Baseline Data Re- 
duction System (Ble cha et al.|2000| ), standard IRAF tasks and cus- 
tom python routines. The sky spectrum was estimated from IFU 
spectra outside an aperture of r ^ 7 spaxels around the object's 
centre and the dedicated sky fibres placed around the IFU. The 
data were corrected for heliocentric velocities and the observations 
from all nights were resampled to a common logarithmic wave- 
length scale and organised into Sd data cubes. For each data cube, 
the centre of the UCD and the effective seeing were determined 
by fitting a model image to the reconstructed image. The model 
image consisted of an ACS high resolution channel (HRC) im- 
age in the F606W filter (shown in the upper left panel of Fig. [T] 
HST programme 101S7, [Evstigneeva et al.|2007] ), convolved with 
a Gaussian PSF and rebinned to the IFU's spatial scale. The five 
observations with the best seeing (0.''50 ^ PSF FWHM ^ Of'67, 
with a mean of 0%0), were then stacked by resampling them to 
a common coordinate grid, at the same time removing the effect 
of atmospheric dispersion { |Filippenko|1982|). The res ampling was 
done using the drizzle algorithm ( [Fruchter & Hook||2002j ), treat- 
ing each spatial slice as an image and masking bad spaxels (cf. the 
weight map in Fig. [T]). These comprised a hot column on the de- 
tector (see the Flames User Manual), as well as spectra that were 
strongly affected by the stray light from neighbouring simultaneous 
calibration lamp spectra. To minimise the loss of spatial resolution 
due to resampling, a moderate amount of super- sampling was ap- 
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plied by using output pixels smaller by a factor of 1.5 on each side. 
The corresponding reconstructed image is shown in Fig.[T| 

The spatial extent of the UCD, with a half-light radius of 
Q'.'95, is only ~2-3 times as large as the spaxel size and the see- 
ing HWHM. In order to preserve the structure of the object, while 
achieving a higher signal to noise ratio in the outer parts of the 
UCD, we used a geometrical scheme for binning the limited num- 
ber of spaxels receiving significant flux from the UCD: The spaxels 
were grouped into radial annuli, which were subdivided into quad- 
rants (shown as red polar grid overlaid on the reconstructed image 
in Fig.[TJ. 



3 KINEMATIC MEASUREMENTS 

The line-of- sight velocity v and velocity dispersion o were recov- 
ered from each spectrum using the penalised pixel-fitting (ppxf) 
code of |Cappellari & Emsellem| ( |2004| ). Two exemplary spectra 
with their best-fitting models are shown on the right of Fig. [T] As 
spectral templates for the fitting we used spectra of field stars and 
members of the NGC 6475 open cluster from the UVES Paranal 
Observatory Project (UVESPOP, [Bagnulo et al.|2003| ) library. The 
library consists of reduced and calibrated high-resolution {R ^ 
80000) stellar spectra spanning a range of spectral- types and metal- 
licities. We used a sub- sample of 46 late-type stars (spectral types 
FO to Ml, luminosity classes II- V), of which the ppxf algorithm 
would select a subset of < 8 significant templates during the fit- 
ting. The spectra were brought to a common radial velocity and 
rebinned to the same logarithmic wavelength- scale as our ARGUS 
data. To degrade the template library to the resolution of our ob- 
servations, we constructed an empirical model of the spectral line- 
spread function (LSF) as a function of wavelength (e.g. |Chilingar-| 
|ian et al.||2007| l, by fitting UVESPOP templates of the according 
spectral type to late- type stars observed with the ARGUS LR04 
setup, and convolved the templates with this LSF. 

The magnitude of the noise in our spectra varied gradu- 
ally with wavelength. Therefore it was modelled by using the 
DER_SNR algorithm ( Stoeh r et al.|200 7 ) to estimate the noise at 
each wavelength and then fitting a linear function to these esti- 
mates. The resulting noise spectrum (shown as green lines in the 
spectra in Fig. [TJ was used as an input for ppxf to ensure the cor- 
rect relative weighting of pixels in the -minimisation, and as a 
basis for generating 100 Monte Carlo realisations of each spectrum, 
which were used to obtain uncertainties on the measured velocity 
and velocity dispersion. 

The fitting was done in the wavelength range from 5156 A to 
5691 A. We masked out spectral regions containing a strong emis- 
sion line doublet in some of the spaxels, most likely originating 
fro m a background galaxy s lightly offset from the UCD's position 
(cf. Evstigneeva et al.|2007| ), the 5577 A 01 sky line, as well as re- 
gions particularly affected by the calibration lamp stray light. For 
consistency the same mask (represented by the grey wavelength in- 
tervals in Fig. [T]) was applied to all spectra, regardless of whether 
they were affected by these blemishes or not. 



4 INTERNAL DYNAMICS 
4.1 Rotation 
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Figure 2. The signature of rotation in the UCD. Left: The top panel shows 
the measured velocity versus the luminosity-weighted azimuth ^ of each 
bin in the outer annulus (r^ I'.' 2) of the spatial binning shown in Fig. ^ 
(squares) and the same binning with the quadrants rotated by 45 degrees 
(circles). While the orientation of the quadrants in Fig.[^is chosen such that 
the maximum rotational signal is captured, the signal essentially disappears 
when the quadrants are rotated by 45 degrees. The blue curve represents 
the best-fitting sinusoid (see Sect. |4.1^ to the velocities measured in both 
orientations of the quadrants (i.e. square and circle symbols). The bottom 
panel shows the same for the inner annuli (r< I'.' 2). Right: The rotation 
curve of the UCD along the direction approximately perpendicular to the 
kinematic axis. 



recover the projected kinematic axis, we initially varied the po- 
sition angle of the quadrants in the spatial binning; once its ap- 
proximate orientation was found, a sinusoid of the form v(r) = 
vq +Asin(0 — PA) was fitted simultaneously to the velocities re- 
covered from that orientation and from the quadrants rotated by 45 
degrees (Fig.[2|. Here vq is the systemic velocity, A is the amplitude 
of the rotation, is the azimuth of each bin, and PA is the position 
angle of the rotational axis. In the outer annulus (r^ X'.'l; upper 
left panel of Fig.|2}, the fit yielded a PA of 12 ± 15 de grees and an 
amplitude of 2.8 ± 0.7 km s~^. Due to the unknown inclination of 
the rotational axis and due to the fact that the maximum of the ro- 
tation curve velocity need not be reached at r~ 1.''5, this represents 
a lower limit on the true rotational velocity. 

4.2 Dispersion profile 

The radial dispersion profile of UCD3 is shown in Fig. [3] The dis- 
persion decreases continuously with increasing radius. There is no 
sign of a disturbed profile, which could have been expected if the 
UCD had undergone a recent interaction, leaving it out of dynami- 
cal equilibrium, a scenario proposed to explain the high dynamical 
masses of some UCDs (Fellhauer & Kroupa 2006). 

Along with the observed velocity dispersion profile. Fig. |3] 
shows a number of models. The models are based on the best-fitting 
two-component luminosity distribution determined by Evstigneeva 



et al.| (|20071 ), a Fornax cluster distance of 18.97 Mpc (Freedman 
et al. [2001] ) and an isotropic velocity distribution. They were con 



The resulting velocity field of the UCD shows the signature of weak 
rotation (see the second panel in the bottom row of Fig. [TJ. To 



structed using the code described in |Hilker et al.|p007| ) and adapted 
to integral-field data. The code takes into account the convolution 
with the seeing PSF and predicts the full line-of- sight velocity dis- 
tribution (LOSVD) in each of our radial annuli. While the mea- 
surement of the observed velocity dispersion assumed a Gaussian 
LOSVD, the actual LOSVD need not be strictly Gaussian. In the 
case of the models including a massive black hole (Sect. |4.3.2[ ), 
where a significant fraction of the observed stars move in the po- 
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agreement with the stellar population M/Ly of 3.7 ±0.2 ( |Chilin-| 
[garian et al.| 20ir). However, regarding the stellar population M/L 
of UCD3, we note that there is a disagreement on the age of UCD3 
in the literature: Lick/IDS absorption indices suggest an age of 
- 2 - 5 Gyr ^Mieske et al.|[2006l [Firth et al.p009 ). This would 
render UCD3 one of the youngest known UCDs and imply a stellar 
M/Ly < 2, which in turn would still require substantial dark mass 
according to our measured dynamical M/L ratio. On the other hand, 
studies based on spectral fitting yield ages of > 12 Gyr ( Chilingar-| 
|ian et aL|[2008l [25TT] ), placing UCD3 within the typical range of 
ages found for most UCDs, which would make our dynamical M/L 
ratio fully consistent with a standard stellar mass function. As a 
possible cause for the discrepant ages, [Firth et al.|p009 ) point out 
that the young age inferred from absorption indices could be an 
artifact introduced by the presence of hot horizontal branch stars. 



Figure 3. The measured velocity dispersions are shown as black squares. 
The radii of the data points correspond to the luminosity-weighted radius 
of each spatial bin. Models based on an isotropic velocity distribution are 
represented by open circles (evaluated in each radial annulus), connected by 
curves to guide the eye: A 'mass follows light' model (blue, solid curve), 
two models including a massive black hole of 5 (red) and 20 per cent (green, 
dashed curves) of the UCD's mass and two models including a dark matter 
component contributing 33 (magenta) and 60 per cent (cyan, dotted curves) 
of the total mass within a three-dimensional radius of 200 pc. 



tential of the central point source, the integrated velocity distribu- 
tion can even become highly non-Gaussian (cf. |Merritt et al.|2009] ). 
Therefore, to obtain the model velocity dispersions that correspond 
to the full modelled LOSVDs produced by the model, we created 
artificial spectra by convolving UVESPOP spectra with the full 
LOSVDs and then measured the dispersion in the same way as for 
the observed spectra. 

Using the fact that the velocity dispersion a scales with the 
mass as (7^ oc M (virial theorem), the best-fitting mass (or equiva- 



lently the best-fitting M/L, for a given observed luminosity L) for 
ea ch model was calc ulated by scaling the model to the data (eq. 5 1 

vne 



Cappellari 



2008 



so that f - = I/(cTm,i-CJobs,i)V(A 



obs.i^ 



minimisea, where CTobs,i ^re the observed velocity dispersions, 
Aobs,i are their associated uncertainties and (7obs,i are the model 
velocity dispersions at the same radii. Confidence regions around 
the best-fitting parameters were estimated from the increase in the 
reduced = k~^x^ compared to its minimum value, where k 
( luiiiibcr of data points — number of model parameters — 1), with 
the number of model parameters being 1 for the constant M/L 
model and 2 for the models including either dark matter or a black 
hole. For one parameter of interest, the 68, 90 and 99 percent confi- 
dence regions are given by an increase Ax^ = ~ Xmin 1^ 2.71 



2007}. For an arbitrary Ax , the confidence 
level can be found by numerical integration of the x^ distribution 
function. The minimum reduced x'^ for the best-fitting model, as 
well as the Ax^ for the other models along with the probabilities to 
find the true model parameter at this Ax^ are given in the legend of 
Fig.|3] 

The simplest possible model, assuming that mass follows light 
(or equivalently, a constant M/L; shown as blue, solid curve) pro- 
vides an excellent fit to the observed dispersion profile, imply- 
ing a UCD mass of M=8.2 ± 0.7 x 10^ M0 and a M/L ratio of 
M/Ly=3.6 ± 0.3. The uncertainty corresponds to the l-cr confi- 
dence interval in the sense of the previous paragraph, obtained by 
varying M/L until x'^ increased by Ax^ = 1 • The M/L is in good 



4.3 Constraints on the non-luminous matter content 

4.3.1 Dark matter 

As a possible explanation for the, on average, high dynamical mass 
to light ratios of UCDs, a substantial dark matter content is dis- 
cussed in the literature (e.g. Hageg an et al.|[2005, ,Baumgardt &| 
|Mieske||2008| ). A potential formation channel would then be the 
tidal stripping of nucleated galaxies, leaving behind a remnant nu- 
cleus with properties comparable to those of a UCD { [Bekki et"aL] 
|2003| ). Simulations show that the resulting UCDs can be dark mat- 
ter dominated, if the central concentration of dark matter was en- 
hanced prior to the stripping through, for example, the adiabatic 
infall of gas ( |Goerdt et al.|2008[[Baumgardt & Mieske|2008| ). 

In order to test this scenario, we included in the mass model 
a dark matter component with a density profile mimicking the dark 
matter distribution predicted for the final evolutionary stage by the 
simulations of |Goerdt et al.| ( [2008| ). Their Fig. 5 shows a final 
dark matter density approximately proportional to r~^-^ for radii 
200 pc< r < 2kpc. We chose as an analytic expression in our model 
a profile of the form p = po(l + r^/r^)~^"^^, which shows the same 
slope at large radii and converges to a central density po at r = 0. 
As a scale radius we adopted = 200 pc being the approximate 
radius below which the |Goerdt et al.| ( [2"008 ') profile starts to flatten. 
We calculated models for a fraction of 20, 33, 60 and 67 per cent 
of dark matter within a three-dimensional radius of 200 pc. We find 
that, while yielding a worse fit than the model without dark mat- 
ter, a dark matter fraction of 33 per cent (shown as magenta dotted 
curve in Fig. [3} is, with Ax^ = 1.3 (again for one parameter of inter- 
est, we marginalise over the M/L), compatible at the l-cr level with 
the observed velocity dispersion profile and yields a stellar M/Ly 
of 3.4. A model with a dark matter fraction of 60 per cent (cyan 
dotted curve) and a best-fitting stellar M/Ly of 2.7, is excluded by 
the data with 98 per cent confidence (Ax^ = 5.5). 



4.3.2 Massive black hole 

Unseen mass in UCDs could also be present in the form of a mas- 
sive black hole. In massive galaxies, nuclear star clusters can co- 
exist with central black holes with comparable masses ( [Graham &| 
|Spitler|2009| ). Therefore, if UCDs represent the remnant nuclei of 
originally much more massive galaxies, at least some of them could 
contain a black hole which contributes a substantial fraction of the 
total mass. Also if UCDs represent compact star clusters around 
recoiling super-massive black holes that were ejected from the cen- 
tres of massive galaxies ( ^Merritt et al.|2009 ), they would contain a 
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black hole with a mass comparable to that of the stellar component. 
Therefore, we explored a family of mass models which, in addi- 
tion to the constant M/L stellar component, included massive black 
holes with 5, 10, 15 and 20 per cent of the UCDs total mass. 

The resulting model dispersion profiles, for a black hole 
mass of 5 (red) and 20 per cent (green) are shown as dashed 
curves in Fig. [3] The models yield best-fitting stellar M/Ly of 2.7 
(MBH/Mtotai=5 per cent) and 1.8 (MBH/Mtotai=20 per cent) and pro- 
vide a worse fit to the data than the model without a black hole, with 
MBH/Mtotal =5 per cent still being consistent with the data at the 1- 
<7 level (A;^^ = 0.8), and the MBn/Mtotai =20 per cent model being 
excluded with 96 per cent confidence (A;^^ =4.3). 



5 CONCLUSIONS 

We have demonstrated that the resolution of the internal kinematics 
of the brightest and most extended UCDs is feasible with seeing- 
limited integral-field spectroscopy, making it a promising tool to 
constrain the formation scenarios and hence the nature of UCDs. 
Our analysis of UCDS in the Fornax cluster reveals weak rotation 
and a dispersion profile in excellent agreement with an isotropic 
velocity distribution and a constant M/L. The mass of the UCD, 
measured under these assumptions, is 8.2 ± 0.7 x 10^ M0 and the 
implied M/L is in good agreement with the most recent determina- 
tion of the UCD's stellar population age and metallicity. The UCD 
does not show any of the attributes, such as tidal disturbance, a 
black hole with a mass comparable to that of the UCD, or a sig- 
nificant dark component, that would strongly suggest that the UCD 
is the remnant nucleus of a larger galaxy or a remnant star cluster 
around a recoiling super-massive black hole. While this result does 
not exclude such scenarios, the internal kinematics of UCDS are 
fully consistent with it being a massive globular cluster. 



ACKNOWLEDGEMENTS 

We thank the anonymous reviewer for valuable comments. MJF 
would like to thank Marina Rejkuba, Eric Emsellem and Harald 
Kuntschner for helpful discussions about this project. HB acknowl- 
edges support from the Australian Research Council through Fu- 
ture Fellowship grant FT0991052. 



REFERENCES 

Bagnulo S., Jehin E., Ledoux C, Cabanac R., Melo C, Gilmozzi 
R., The ESO Paranal Science Operations Team 200S, ESO Mes- 
senger, 114, 10 
Baumgardt H., Mieske S., 2008, MNRAS, S91, 942 
Bekki K., Couch W. J., Drinkwater M. J., Shioya Y., 200S, MN- 
RAS, S44, 399 

Blecha A., Cayatte V., North R, Royer F, Simond G., 2000, in 
M. lye & A. F Moorwood ed., SPIE Conf. Ser. Vol. 4008, Data- 
reduction software for GIRAFFE, the VLT medium-resolution 
multi-object fiber-fed spectrograph, p. 467 
Cappellari M., 2008, MNRAS, S90, 71 
Cappellari M., Emsellem E., 2004, PASP, 116, 1S8 
Chilingarian I. V., Cayatte V., Bergond G., 2008, MNRAS, S90, 
906 

Chilingarian I. V., Mieske S., Hilker M., Infante L., 2011, MN- 
RAS, pp 66-+ 



Chilingarian I. V., Prugniel P., Sil'Chenko O. K., Afanasiev V. L., 

2007, MNRAS, S76, lOSS 
Da Rocha C, Mieske S., Georgiev I. Y., Hilker M., Ziegler B. L., 

Mendes de Oliveira C, 2011, A&A, 525, A86-h 
Dabringhausen J., Fellhauer M., Kroupa P, 2010, MNRAS, 403, 

1054 

Dabringhausen J., Hilker M., Kroupa P, 2008, MNRAS, 386, 864 
Dabringhausen J., Kroupa P, Baumgardt H., 2009, MNRAS, 394, 
1529 

De Propris R., Phillipps S., Drinkwater M. J., Gregg M. D., Jones 
J. B., Evstigneeva E., Bekki K., 2005, ApJ, 623, L105 

Drinkwater M. J., Jones J. B., Gregg M. D., Phillipps S., 2000, 
PAS A, 17, 227 

Eisenhauer F. et al., 2003, in M. lye & A. F. M. Moorwood 
ed., SPIE Conf. Ser. Vol. 4841, SINFONI - Integral field spec- 
troscopy at 50 milli-arcsecond resolution with the ESO VLT, p. 
1548 

Evstigneeva E. A. et al., 2008, AJ, 136, 461 

Evstigneeva E. A., Gregg M. D., Drinkwater M. J., Hilker M., 

2007, AJ, 133, 1722 
Fellhauer M., Kroupa P, 2002, MNRAS, 330, 642 
Fellhauer M., Kroupa P, 2006, MNRAS, 367, 1577 
Filippenko A. V, 1982, PASP, 94, 715 

Firth P, Evstigneeva E. A., Drinkwater M. J., 2009, MNRAS, 394, 
1801 

Freedman W. L. et al., 2001, ApJ, 553, 47 

Fruchter A. S., Hook R. N., 2002, PASP, 114, 144 

Goerdt T., Moore B., Kazantzidis S., Kaufmann T., Maccio A. V, 

Stadel J., 2008, MNRAS, 385, 2136 
Graham A. W, Spitler L. R., 2009, MNRAS, 397, 2148 
Ha§egan M. et al., 2005, ApJ, 627, 203 

Hau G. K. T., Spitler L. R., Forbes D. A., Proctor R. N., Strader 
J., Mendel J. T., Brodie J. P, Harris W. E., 2009, MNRAS, 394, 
L97 

Hilker M., Baumgardt H., Infante L., Drinkwater M., Evstigneeva 

E., Gregg M., 2007, A&A, 463, 119 
Hilker M., Infante L., Vieira G., Kissler-Patig M., Richtler T., 

1999, A&AS, 134, 75 
Kaufer A., Pasquini L., Castillo R., Schmutzer R., Smoker J., 

2003, ESO Messenger, 113, 15 
Kruijssen J. M. D., Mieske S., 2009, A&A, 500, 785 
McGregor P J. et al., 2003, in M. lye & A. F. M. Moorwood ed., 

SPIE Conf. Ser. Vol. 4841, Gemini near-infrared integral field 

spectrograph (NIFS), p. 1581 
Merritt D., Schnittman J. D., Komossa S., 2009, ApJ, 699, 1690 
Mieske S., Hilker M., Infante L., 2004, A&A, 418, 445 
Mieske S., Hilker M., Infante L., Jordan A., 2006, AJ, 131, 2442 
Mieske S. et al., 2008, A&A, 487, 921 
Mieske S., Kroupa P, 2008, ApJ, 677, 276 
Pasquini L. et al., 2002, ESO Messenger, 110, 1 
Paudel S., Lisker T., Janz J., 2010, ApJ, 724, L64 
Press W. H., ed. 2007, Numerical recipes, 3rd edn.Cambridge 

University Press, Cambridge 
Rejkuba M., Dubath P, Minniti D., Meylan G., 2007, A&A, 469, 

147 

Stoehr F et al., 2007, STECF, 42, 4 

Taylor M. A., Puzia T. H., Harris G. L., Harris W. E., Kissler-Patig 
M., Hilker M., 2010, ApJ, 712, 1191 

This paper has been typeset from a TgK/ LTgX file prepared by the 
author. 



© 2011 RAS, MNRAS 000,[l]|5] 



